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Abstract 

Column-convex polyominoes were introduced in 1950's by Temperley, 
a mathematical physicist working on "lattice gases". By now, column- 
convex polyominoes are a popular and well-understood model. There 
exist several generalizations of column-convex polyominoes; an example 
is a model called multi- directed animals. In this paper, we introduce a 
new sequence of supersets of column-convex polyominoes. Our model 
(we call it level m column- subconv ex polyominoes) is defined in a sim- 
ple way. We focus on the case when cells are hexagons and we compute 
the area generating functions for the levels one and two. Both of those 
generating functions are complicated g-series, whereas the area generat- 
ing function of column-convex polyominoes is a rational function. The 
growth constants of level one and level two column-subconvex polyomi- 
noes are 4.319139 and 4.509480, respectively. For comparison, the growth 
constants of column-convex polyominoes, multi-directed animals and all 
polyominoes are 3.863131, 4.587894 and 5.183148, respectively. 
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1 Introduction 

The enumeration of polyominoes is a topic of great interest to chemists, physi- 
cists and combinatorialists alike. In chemical terms, any polyomino (with hexag- 
onal cells) is a possible benzenoid hydrocarbon. In physics, determining the 
number of n-celled polyominoes is related to the study of two-dimensional per- 
colation phenomena. In combinatorics, polyominoes are of interest in their own 
right because several polyomino models have good-looking exact solutions. 

Known results for polyominoes include the fact that the number of n-celled 
polyominoes grows exponentially. More precisely, if a n denotes the number of 
n-celled polyominoes, then 

n . l/n 1/n 

• lim„_ ) . 00 a n ' = t = sup„> x a n ' 

• lim^oo a n+ i/a n = r. 

The first result follows from standard concatenation arguments, see e.g. [T2HT3"] . 
while the second result, due to Madras [TO] relies on a pattern theorem for lattice 
animals. These results are quite general, and apply mutatis mutandis to the new 
polyomino models we consider here. 

One can also obtain rigorous bounds on the growth constant r. For example, 
for hexagonal polyominoes we have [12] 

4.8049 < r w 5.183148 < 5.9047. 

A lower bound is immediately obtainable from the first itemised equation above, 
and it can be improved with rather more work. The upper bound is obtained by 
a method due to Klarncr and Rivest [9] , which relies on mapping each polyomino 
onto a tree on the dual lattice, and relaxing the rules for tree construction so 
that over-counting results. 

One very popular polyomino model is that of column- convex polyominoes. 
Column-convex polyominoes with hexagonal cells have a rational area generat- 
ing function. That generating function was found by Klarner in 1967 [5]. The 
growth constant of hexagonal-celled column-convex polyominoes is 3.863131. 
(By the growth constant we mean the limit lim^oo ^/a~^, where a n denotes the 
number of n-celled elements in a given set of polyominoes.) 

In a previous paper [4], one of us (Feretic) began to search for polyomino 
models which are more general than column-convex polyominoes, but still have 
reasonably simple area generating functions. In [4], Feretic introduced level 
m cheesy polyominoes (m = 1, 2, 3, . . .), and here we shall introduce another 



2 



sequence of models, which we call level m column- subconv ex polyominoes (to = 
1, 2, 3....B 

At every level, cheesy polyominoes have a rational area generating function, 
whereas column-subconvex polyominoes have an area generating function which 
is unlikely to be algebraic, and indeed, unlikely to be differentiably finite [5]. 
Further, at any given level, cheesy polyominoes are an exponentially small subset 
of column-subconvex polyominoes. The latter set of polyominoes has a greater 
growth constant than the former set. For example, the growth constant of 
level one cheesy polyominoes is 4.114908, while the growth constant of level 
one column-subconvex polyominoes is 4.319139. In addition, if we reflect a 
column-subconvex polyomino about a vertical axis, we get a column-subconvex 
polyomino again. This kind of invariance under reflection is enjoyed by column- 
convex polyominoes, but not by cheesy polyominoes. Admittedly, counting level 
to column-subconvex polyominoes requires more effort than counting level to 
cheesy polyominoes. Anyway, at level one, column-subconvex polyominoes are 
not very hard to count. Just as with cheesy polyominoes, as the level increases, 
the computations quickly increase in size. 

In this paper, the level one column-subconvex model is solved in full detail. 
We also solved the level two column-subconvex model. To see the level two 
result (stated with no proof), the reader may visit the web page [2]. The said 
result involves too many auxiliary expressions to be stated in this paper. (To 
be specific, there are 33 auxiliary expressions, of which 25 are polynomials; the 
degree of those polynomials is between 20 and 23.) 

Our computations are done by using Bousquet-Melou's [1] and Svrtan's [6] 
"turbo" version of the Temperley method pj] . 

If the reader would like to have more information on the history of polyomino 
enumeration, or on the role which polyominoes play in physics and chemistry, 
then he/she may refer to Bousquet-Melou's habilitation thesis [2j, or to the book 

2 Definitions and conventions 

There are three regular tilings of the Euclidean plane, namely the triangular 
tiling, the square tiling, and the hexagonal tiling. We adopt the convention 
that every square tile or hexagonal tile has two horizontal edges. In a regular 
tiling, a tile is often referred to as a cell. A plane figure P is a polyomino if P is 
a union of finitely many cells and the interior of P is connected. See Figure 1. 
Observe that, if a union of hexagonal cells is connected, then it possesses a 
connected interior as well. 

Let P and Q be two polyominoes. We consider P and Q to be equal if and 
only if there exists a translation / such that f(P) = Q- 

If a polyomino P is made up of n cells, we say that the area of P is n. 

1 Bousquet-Melou and Rechnitzer's multi- directed animals [3] are also a superset of column- 
convex polyominoes with hexagonal cells. 
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Figure 1: A hexagonal-celled polyomino. 



Let R be a set of polyominoes. By the area generating function of R we 
mean the formal sum 



From now on, we concentrate on the hexagonal tiling. When we write "a 
polyomino" , we actually mean "a hexagonal-celled polyomino" . 

Given a polyomino P, it is useful to partition the cells of P according to their 
horizontal projection. Each block of that partition is a column of P. Note that 
a column of a polyomino is not necessarily a connected set. An example of this 
is the highlighted column in Figure 1. On the other hand, it may happen that 
every column of a polyomino P is a connected set. In this case, the polyomino 
P is a column- convex polyomino. See Figure 2. 

Let a be a column of a polyomino P. By the height of a we mean the number 
of those cells which make up a plus the number of those (zero or more) cells 
which make up the gaps of a. For example, in Figure 1, the highlighted column 
has height 7, and the next column to the left has height 4. 

A finite union of cells P is a level m column- subconv ex polyomino if the 
following holds: 

• P is a polyomino, 

• every column of P has at most two connected components, 

• if a column of P has two connected components, then the gap between 
the components consists of at most m cells. 




f area of P 



Pen, 
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Figure 2: A column-convex polyomino. 




The gray cell 
is the pivot cell. 



Figure 4: The pivot cell. 



See Figure 3. 

Let S denote the set of all level one column-subconvex polyominoes. 

Let P be an element of S and let P have at least two columns. Then we 
define the pivot cell of P to be the lower right neighbour of the lowest cell of 
the second last column of P. See Figure 4. Observe that the pivot cell of P is 
not necessarily contained in P. 

When we build a column-convex polyomino from left to right, adding one 
column at a time, every intermediate figure is a column-convex polyomino itself. 
However, when we build a column-subconvex polyomino, this is no longer the 
case. A "left factor" of an element of S need not itself be a polyomino, and 
therefore need not be an element of S. 

We say that a figure P is an incomplete level one column-subconvex poly- 
omino if P itself is not an element of S, but P is a "left factor" of an element 
of S. Notice that, if P is an incomplete level one column-subconvex polyomino, 
then the last {i.e., the rightmost) column of P necessarily has a hole. 

Let T denote the set of all incomplete level one column-subconvex polyomi- 
noes. 

Let P be an element of SUT and let P have at least two columns. Then we 
define the body of P to be all of P, except the rightmost column of P. 

Let P be an element of T and let P have at least two columns. We define 
the lower pivot cell of P to be the lower right neighbour of the lowest cell of the 
second last column of P. In addition, we define the upper pivot cell of P to be 
the upper right neighbour of the highest cell of the second last column of P. 

3 Notations for generating functions. Partitions 
of the sets 
S and T 

We shall deal with the following generating functions: 
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J^^q f'j ^ ^ qdrea of P _ ^the height of the last column of P 

Pes 



A x =A(q,\), B 1 = —(q,l), 



C(q,u,v) = ^l^^ 



the height of the upper 
P ^ component of the last column of P 



pet 



the height of the lower 
component of the last column of P 

•V 



D(u) = C(q,u,l), E(v) = C(q,l,v), d=C(g,l,l). 

Functional equations for the generating functions will be obtained by the 
"divide and conquer" strategy. Namely, now we are going to partition the sets 
S and T. 

Let S a be the set of level one column-subconvex polyominoes which have 
only one column. 
Let 



Sp = {P € S \ S a : the body of P lies in S, the last column of P 
has no hole, and the pivot cell of P is contained in P}, 

S 7 = {P € S \ S a : the body of P lies in S*, the last column of P 
has no hole, and the pivot cell of P is not contained in P}, 

Ss = {P € S \ S a : the body of P lies in S*, and the last column of P 
has a hole}, 

S f _ = {P G S \ S a : the body of P lies in T, and the last column of P 
has no hole} and 
= {P G S \ S a : the body of P lies in T, and the last column of P 
has a hole}. 

The sets S ai Sp, S 7 , S$, S e and 5^ form a partition of 5*. We write A a , Ap, 
Ay, A , A e and Aq for the parts of the series A that come from the sets S a , Sp, 
Ss, S e and S^, respectively. 

We proceed to the set T. We write T a for the set of incomplete level one 
column-subconvex polyominoes which have only one column. Let P € T\T a . 
If the body of P lies in S, then the said body is in contact with just one of the 
two connected components of P's last column. The non-contacting component 
of the last column is located either wholly above or wholly below the second 
last column of P. Let 
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Tp = {P £T \ T a : the body of P lies in S, and the hole of the last 

column of P coincides either with the lower pivot cell of P 

or with the upper pivot cell of P} and 
T 7 = {P £T\T a : the body of P lies in S, and the hole of the last 

column of P lies either below the lower pivot cell of P 

or above the upper pivot cell of P}. 

Let us move on to the case when the body of P £ T \ T a lies in T. Then 
the second last column of P has two connected components. It is easy to see 
that each of those two components must be in contact with the last column 
of P. (This docs not mean that each of the two connected components of 
the last column of P must be in contact with the second last column of P.) 
Now, it may or may not happen that one connected component of P's last 
column is in contact with both connected components of P's second last column. 
Accordingly, we define the following two sets: 



T s = {P £T\T a : the body of P lies in T, and the hole of the last 

column of P touches the hole of the second last column of P} and 

T e = {P £T\T a : the body of P lies in T, and the hole of the last 

column of P docs not touch the hole of the second last column of P}. 

The sets T a , Tp, T 7 , Ts and T e form a partition of T. We write C Q , Cp, C 7 , 
Cs and C e for the parts of the series C that come from the sets T ai Tp, T 7 , Ts 
and T e , respectively. 



4 Setting up the functional equations for A, A\ 
and B 



To begin with, it is clear that 

A a =qt + {qt) 2 + (qt) 3 + . 



qt 



(1) 



1 - qt 

If a polyomino P lies in Sp, then the last column of P is made up of the pivot 
cell, of i £ {0, 1, 2, 3, . . .} cells lying below the pivot cell, and of j £ {0, 1, 2, 3, . . .} 
cells lying above the pivot cell. See Figure 5. Hence, 



Ap = A 1 -qt- 



3=0 



qt 



(1 - q ty 



(2) 



Consider the following situation. A polyomino P £ S ends with a column 
I. We are creating a new column to the right of /, and the result should be an 
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Figure 5: The last two columns of two elements of Sp. 
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Figure 6: The last two columns of two elements of S 7 . 



element of 5 7 . Then, whether or not the column I has a hole, we can put the 
lowest cell of the new column in exactly m places, where m is the height of /. 
See Figure 6. Hence 

Ar^-B, (3) 

Let us proceed to another situation. A polyomino P e S ends with a column 
J. We are creating a new column to the right of J, and the result should be an 
element of Sg. Then, whether or not the column J has a hole, we can put the 
hole of the new column in exactly n — 1 places, where n is the height of J. See 
Figure 7. The new column is made up of i e {1, 2, 3, . . . } cells lying below the 
hole, of a hole of height one, and of j e {1, 2, 3, . . . } cells lying above the hole. 
Altogether, 

*-iV-i^-<*-*>-(i^-<*-*>- (4 ' 

Now, let P be an element of S e . By the definition oi S e , P is a polyomino 
with a one-part last column, but the body of P is not a polyomino. Hence, 
in the second last column of P there is a hole, and in the last column of P 
there are two cells with which the hole is filled. In addition to this two-celled 
"cork", the last column contains i € {0, 1, 2, . . . } cells lying below the "cork" 
and j 6 {0, 1, 2, ... } cells lying above the "cork". See Figure 8. Hence 
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Figure 7: The last two columns of two elements of S$. 
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Figure 8: The last two columns of an element of S e . 



A = 



l 




l 



■c x = 




(5) 



l-qt 



l-qt 



(i - q ty 



If P is an element of Sq, then P is a polyomino with a two-part last column, 
while the body of P is not a polyomino. Once again, in the second last column 
of P there is a hole, and in the last column of P there are two cells with which 
the hole is filled. Let the lower component of the second last column consist of 
i cells, and let the upper component of the second last column consist of j cells. 
Now, if the two-celled cork belongs to the upper component of the last column, 
then it is impossible that i — 1 or more cells of the last column lie between the 
cork and the hole of the last column. Namely, if i — 1 or more cells were so 
situated, then the lower component of the last column would not be connected 
with the rest of P, and P would not be a polyomino. See Figure 9. 

For a similar reason, if the two-celled cork belongs to the lower component of 
the last column, then it is impossible that j — 1 or more cells of the last column 
lie between the cork and the hole of the last column. 

These remarks lead us to the following expression: 
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Figure 9: The last two columns of an element of S$. 
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n 3 f 4 3 f 4 



(1-qtf {1-qtf 
aH 4 o 3 t 4 



When an element of T is reflected about a horizontal axis, the area is 
preserved, whereas the height of the upper component of the last column be- 
comes the height of the lower component of the last column. This means that 
D(u) = E(u). We now have 

Since A = A a + Ap + A 7 + A s + A e + A$, equations (4.1)-(4.6) imply that 
qt qt qt q 2 t 3 

A —t + Ti^- Al + —t- Bl + T^—W- {Bl - Al) 

qH 2 2qH 4 2fl 



(1 - qt) 2 (1 - qtf (1 - qtf 

Setting t = 1, from equation (4.7) we get 

Differentiating equation (4.7) with respect to t and then setting t = 1, we 

get 

q q + q 2 _q 

+ (i - ^• (i?1 " yll) + (i _ ^' Cl + 7r^F' Cl 

-(T^-^)-(T^-^)- («) 

5 Setting up the functional equations for C, D 
and Ci 

Now we turn to incomplete level one column-subconvex polyominoes. We have 
already observed that an incomplete level one column-subconvex polyomino 
always ends with a holed column. 
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Figure 10: The last two columns of two elements of Tp. 



The set T a contains every two-part column (with one-celled hole) having 
i e {1, 2, 3, . . . } cells below the hole and j e {1, 2, 3, ... } cells above the hole. 
Thus, 

c a = qv ■ qu = q2uv (10) 

a 1 — qv 1 — qu (1 — qu)(l — qv) ' 

If P e Tp, then the body of P lies in S. The hole of the last column has two 
possibilities: to coincide with the lower pivot cell of P or to coincide with the 
upper pivot cell of P. Anyhow, the last column is made upofiG {1, 2, 3,...} 
cells lying below the hole and j e {1, 2, 3, ... } cells lying above the hole. See 
Figure 10. Therefore, 

C/3 = 2 .^.^. Al = ? ^ Al . (n) 

1 — qv 1 — qu (l — qu)(l — qv) 

Now let P £ T 7 . The body of P again lies in S. If the hole of the last 
column lies below the lower pivot cell of P, then the last column of P is made 
up of: 
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• i G {1, 2, 3, ... } cells lying below the hole, 

• j S {0, 1, 2, ... } cells lying above the hole and below the lower pivot cell, 

• the lower pivot cell, and 

• k e {0, 1, 2, . . . } cells lying above the lower pivot cell. 

If the hole of the last column lies above the upper pivot cell of P, then the 
last column of P is made up of: 



i e {1, 2, 3, . . . } cells lying above the hole, 

j G {0, 1, 2, ... } cells lying below the hole and above the upper pivot cell, 
• the upper pivot cell, and 

k e {0, 1, 2, . . . } cells lying below the upper pivot cell. 

See Figure 11. Altogether, 

q 2 uv q 2 uv 
° 7 ~ (1 - gu) 2 (l - qv) 1 + (1 - qu)(l - qvf ' ^ [ > 

If P £T$, then the body of P lies in T. The second last and last columns of P 
both have a hole. The hole of the last column is either the lower right neighbour 
or the upper right neighbour of the hole of the second last column. In the last 
column, there are i 6 {1, 2, 3,...} cells below the hole and j £ {1, 2, 3, ... } 
cells above the hole. See Figure 12. Hence, 

c i = 2 .JL.JL. Cl = - 2 ^ . Cl . (13) 

1 — qv 1 — qu {l — qu)(l — qv) 

Let P &T e . Once again, the second last and last columns of P both have a 
hole. However, to the right of the hole of the second last column, there are two 
cells which both belong to P. If this two-celled cork is contained in the lower 
component of the last column, and if the upper component of the second last 
column consists of j cells, then it is necessary that at least j — 1 cells of the last 
column lie above the cork and below the hole of the last column. Otherwise the 
upper component of the last column would be connected with the rest of P, and 
P would be a polyomino. (That cannot happen because P is an element of the 
set T, and the elements of T are not polyominoes.) See Figure 13. 

Similarly, if the two-celled cork is contained in the upper component of the 
last column, and if the lower component of the second last column consists of 
i cells, then it is necessary that at least i — 1 cells of the last column lie below 
the cork and above the hole of the last column. If it were not so, then P would 
be a polyomino. 

These remarks made, we conclude that 
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Figure 11: The last two columns of two elements of T 7 . 
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Figure 12: The last two columns of two elements of T$. 
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Figure 13: The last two columns of the elements of T e . 
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(1 — qu)(l — qv) 2 (1 — qu) 2 (l — qv) y 

D (<P>) + T, \rr, s-D(qu). (14) 



(1 - qu){\ - qv) 2 ^' (l-qu) 2 (l-qv) 
Since C = C a + Cp + C 7 + Cs + C £ , equations (5.1)-(5.5) imply that 



g to 2q uv ^ q uv ^ 

(1 — qu)(l — qv) (l — qu)(l — qv) 1 (1 — qu) 2 (l — qv) 1 

H g 2 ^ A 2q 2 uv c 

(1 — qu)(l — qv) 2 (l — qu)(l — qv) 

9 9 

q uv q uv 

+ (l-qu)(l-qv) 2 ■ D{qV) + (l- q u) 2 (l-qv) ' (15) 



Setting u = 1, from equation (5.6) we get 



2 2 

(U) = (l- q )(l-qu) + (l-q)(l-qu) 2 ' 1 

I (3 ~ 2g)g 2 u 2g 2 M 
+ (1- 1 + (1 - q)(l - qu) 1 

+ ^-q%-quY D ^ ^-q^-quf - D ^ <»> 

We have D(l) = C(q, 1,1) = Ci. So, when we set u = 1, equation (5.7) 
turns into 



6 Solving the functional equations 

For convenience, we firstly define an extra series F. The definition is 

F=l + ^ r ^--A 1 +2C 1 + T ^-D(q). (18) 
Now equation (5.7) can be written as 



_., . q 2 u q 2 u q 2 u 

D(U) = (1 - 9 )d - .u) 2 ' Al + (1 - q)(l - qu) (I - q)(l - qu) 2 ' 

(19) 
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The next step of the upgraded Temperley method is to iteratively remove 
the D(qu) term from the right-hand side of equation (6.2). Namely, substituting 
qu for u in equation (6.2) gives 



D(qu) 



q 3 u 



(l-#-9 2 «) : 



q 3 u 



■F 



q 3 u 



{l-q){l-q 2 u)- 



■D(q 2 u). 



Replacing D{qu) of equation (6.2) by this latter expression, followed by a 
bit of rearranging, results in 



D(u) 



9 

q U 



q 2+3 u 2 



{l-q){l-qu) 2 {l-q) 2 {l-qu) 2 {l^q 2 uf 



9 

q z u 



q 2+3 u 2 



(l-q)(l-qu) {I- q) 2 {\- qu) 2 {l- q 2 u) 

D{q 2 u). 



■A 1 
■ F 



q 2+3 u 2 



{l-q) 2 {l-qu) 2 {l-q 2 uf 



After the next iteration, in each of the square brackets there is a sum of three 
terms, and the argument of the hnal D is q 3 u instead of q 2 u. After infinitely 
many iterations, we have 



D(u) 



E 



q 2 u 



A, 



q 2 u 



E 

<=i (i~ q y- rpfc=i(i-? fc «) .(i-g*«) 



F. (20) 



The right-hand side of equation (6.3) involves no D because lim n ^ 00 D{q n u) — 
0. The reason why this limit is zero is that the lowest power of q occurring in 
D(q n u) is n + 2. 

Setting u = q, from equation (6.3) we get 



w = E rV : 

a-?) 1 - nu(i-9 fe+i ) 



Ax 



E 



1 a-?) 1 - ni=i(i-? 



(i-<f+i; 



F (21) 



"Logarithmically" differentiating! equation (6.3) with respect to u and then 
setting u = q, we obtain 



2 By logarithmic differentiation we mean the use of the formula ip' = ip ■ [ln(<p)]' . 
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d'm = E 

i=i(i-g)*- nu(i 



E 



,fc+i s 



9 fel" 



=1 (i - «)* ■ nl=i(i-9 fe+l ) 



+ 



1 - 



(22) 



Equations (4.8), (4.9), (5.8), (6.1), (6.4) and (6.5) make up a system of 
six linear equations in six unknowns: A\, B\, Ci, D(q), D'(q) and F. That 
linear system was readily solved by the computer algebra package Maple. Of 
course, the most interesting component of the solution is A\, the area generating 
function for level one column-subconvex polyominoes. We state the formula for 
Ai as a theorem. 

Theorem 1 The area generating function for level one column-subconvex poly- 
ominoes is given by 



A l = 



ELi num n 

En=l den « 



where 



num x = q-8q 2 + 2Aq 3 -i2q 4 + \7q b + Aq e -8q 7 + 2q 8 , 

num 2 = {-q + hq 2 - 13g 3 + 23g 4 - 22 9 5 + 12g 6 - 2q 7 ) ■ 0, 

num 3 = (-2q 4 + 8q 5 - I2q 6 + 8q 7 - 2q s ) ■ 5, 

dem = l-llg + 46g 2 - 93q 3 + 88q 4 -27q 5 - 2Aq e + 19q 7 - 3<7 8 , 

den 2 = {2q 2 - 8q 3 + 8q 4 - Aq 5 - 6q 6 + Aq 7 ) ■ a, 

den 3 = {-1 + lOq - 3Aq 2 + 67 q 3 -8lq 4 + 5Aq 5 -16q 6 + q 7 ) ■ 0, 

den A = (2<7 4 - 8q 5 + 8q 6 - 2q 8 ) ■ 7, 

den 5 = (6<7 4 - 22q 5 + 34g 6 - 22q 7 + Aq 8 ) ■ 5, 

den 6 = {2q 4 - 6q 5 + lOq 6 - 6q 7 ) ■ (aS - £7), 



a 



2 ' 
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2 (1-^+1) ' 

2 ' 

From the formula just stated, one easily finds that 

A 1 = g + 3<7 2 + n<? 3 + 44g 4 + 184,? 5 + 786q 6 + 339f + 14683g 8 
+63619g 9 + 275506g 10 + f 192134,? 11 + 51 54794 g 12 + . . . . 

We expanded A\ in a Taylor series to 250 terms, and analysed the series by 
the method of differential approximants [7] using second-order approximants, 
that is to say, approximants given by solutions of inhomogeneous second degree 
ordinary differential equations. From this analysis, we found that the dominant 
singularity of A\ is a simple pole, located at q — q c — 0.2315276132 . (Note 
that we only needed some 20 series terms to establish this — the additional terms 
merely provided higher accuracy and confirmation of our initial analysis). We 
refined this estimate by using Maple to locate the position of the denominator 
zero. That is, by expanding the series Y^n=i den n to more and more terms, 
more and more accurate numerical solutions of Y] n _ t den n — were obtained. 
In this way, we found q c to be 0.231527613159. We could obtain much higher 
accuracy if necessary. It is likely that this number is algebraic (such is usually 
the case with exact solutions), but we have been unable to conjecture its exact 
value. Note also that the numerator is positive for < q < 0.8, so there is no 
possibility that this denominator zero cancels with the numerator. 

The growth constant is the reciprocal of the dominant singularity. Thus, the 
growth constant of level one column-subconvex polyominoes is about 4.319139. 
For comparison, the growth constant of column-convex polyominoes is 3.863131, 
the growth constant of level one cheesy polyominoes is 4.114908, and the growth 
constant of all polyominoes is 5.183148. (The latter two growth constants were 
found in [4] and [12], respectively.) From the result that r = sup n>1 al/ n , we also 
have the quite good lower bound (based on 250 terms) t > 4.283006. We can also 
calculate the amplitude, so writing the generating function as Ai — ^ n a n q n , 
then a n ~ c\ ■ r™, we can estimate the amplitude c\ from the sequence of 
quotients a n /r n . In this way we estimate c\ = 0.1224281005. 
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7 Level two column-sub convex polyominoes 



In the just-finished enumeration of level one column-subconvex polyominoes, we 
considered altogether 11 cases. Namely, we partitioned the set S into 6 subsets 
and the set T into 5 subsets; 6 + 5 equals 11. 

Let us write A 2 to denote the area generating function for level two column- 
subconvex polyominoes. We found a formula for A 2 as well, but that goal was 
achieved through considering as much as 37 cases. We had to struggle against 
a number of complicated expressions, and in the end we had to solve a system 
of 16 linear equations in 16 unknowns. (For comparison, the computation of 
A\ was completed by solving a system of 6 linear equations in 6 unknowns.) 
Consequently, the formula for A 2 is much bulkier than the formula for A\. To 
save this journal's space, we have chosen to state the formula for A 2 in the 
electronic form only [14] . 

However, the Taylor series expansion of A 2 is 

A 2 = g + 3g 2 + ll< ? 3 + 44g 4 + 186<? 5 + 812< ? 6 + 3614g 7 + 16254«? 8 
+ 73464g 9 + 332603g 10 + 1505877<j n + 6813301(7 12 + . . . 

and the critical point of A2 is at q = q c = 0.221755050048. This was obtained in 
the same way as described above for level one column-subconvex polyominoes, 
but based on a series of 153 terms. Thus the growth constant of level two 
column-subconvex polyominoes is about 4.509480. For comparison, the growth 
constant of level two cheesy polyominoes is 4.231836 4J. As above, we can 
also give the rigorous bound r > 4.441222. We can also estimate the amplitude 
C2 = 0.0969488405, so that the coefficient of the n th term of the generating 
function A 2 = J2 n a ™ 1 ?™ is a « ~ c 2 ' T ™- 

As stated in [12j . the area generating function for all polyominoes is 

q + iq 2 + lie? 3 + 44 g 4 + 186g 5 + 814g 6 + 3652g 7 + 16689g 8 
+ 77359<? 9 + 362671<7 10 + 1716033g n + 8182213g 12 + . . . . 

Indeed, a quick drawing confirms that a polyomino must have at least 5 
(resp. 6) cells in order not to be a level one (resp. two) column-subconvex 
polyomino. See Figure 14. 

We have not tried to enumerate level three column-subconvex polyominoes. 
Our non-rigorous estimate is that, in order to enumerate this latter model by 
area, one would have to consider at least 80 cases. 

8 Conclusion 

We have defined a class of polyominoes that interpolates between column-convex 
polyominoes and all polyominoes. The former have been solved, while the latter 
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Figure 14: (a) The two smallest instances of a polyomino which is not a level one 
column-subconvex polyomino. (b) The two smallest instances of a polyomino 
which is not a level two column-subconvex polyomino. 



remain unsolved. For now, our interpolating class (we call it level m column- 
subconvex polyominoes) remains solved up to a certain point. Namely, we have 
solved the cases m = 1 and m = 2. Column-convex polyominoes correspond 
to the case m = 0. In both cases m = 1, 2, the generating function has a 
simple pole singularity, located at q = q c = 0.2315 . . . and 0.2217 . . . respectively. 
For all polyominoes, the corresponding singularity is at q — q c (polyomino) = 
0.192932... , and the singularity is of the form const. ■ |log(g c — q)\, rather 
than a simple pole [12j . For all finite values of m we expect the generating 
function of level m column-subconvex polyominoes to have a simple pole, while 
the singularity position is expected to be a monotone decreasing function of m, 
with a limiting value as m tends to infinity of q* > q c (polyomino). We have also 
given the rigorous lower bounds r > 4.283006 and r > 4.441222 for the growth 
constants of level 1 and level 2 column-subconvex polyominoes respectively. 

References 

[1] M. Bousquet-Melou, A method for the enumeration of various classes of 
column-convex polygons, Discrete Math. 154 (1996), 1-25. 

[2] M. Bousquet-Melou, Rapport Scientifique pour obtenir l'habilitation a 
diriger des recherches, Report No. 1154-96, LaBRI, Universite Bordeaux I, 
1996. 

[3] M. Bousquet-Melou and A. Rechnitzer, Lattice animals and heaps of 
dimers, Discrete Math. 258 (2002), 235-274. 

[4] S. Feretic, Polyominoes with nearly convex columns: A semidirected model. 
Math. Commun., submitted. 

[5] S. Feretic and A.J. Guttmann, Two generalizations of column-convex poly- 
gons, preprint. 



25 



[6] S. Feretic and D. Svrtan, On the number of column- convex polyominoes 
with given perimeter and number of columns, in: Proc. of the Fifth FPSAC 
Conference (eds. A. Barlotti, M. Delest and R. Pinzani), Firenze, 1993, 
201-214. 

[7] A. J. Guttmann, Analysis of coefficients, in: Phase Transitions and Critical 
Phenomena, Vol. 13 (C. Domb and J. L. Lebowitz, eds.), Academic Press, 
New York, 1989, 1-234. 

[8] D. A. Klarner, Cell growth problems, Canad. J. Math. 19 (1967), 851-863. 

[9] D. A. Klarner and R. L. Rivest, A procedure for improving the upper bound 
for the number of n-ominoes, Canad. J. Math. 25 (1973), 585-602. 

[10] N. Madras, A pattern theorem for lattice clusters, Ann. Comb. 3 (1999), 
357-384. 

[11] H. N. V. Temperley, Combinatorial problems suggested by the statistical 
mechanics of domains and of rubber-like molecules, Phys. Rev. 103 (1956), 
1-16. 

[12] M. Voge and A. J. Guttmann, On the number of hexagonal polyominoes, 
Theoret. Comput. Sci. 307 (2003), 433-453. 

[13] Polygons, Polyominoes and Polycubes (ed. A. J. Guttmann), Lecture Notes 
in Phys., Vol. 775, Springer, Berlin, 2009. 

[14] \protect\vrule widthOpt\protect\href {http : //www.gradri .hr/adminmax/f iles/staf f /A_2 . 

(This file is a Maple 9.5 worksheet. The file can also be obtained from 
Svjetlan Feretic via e-mail.) 



26 



